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Abstract 


In  the  analysis  of  current -voltage  characteristics  of  a  probe  in  a  plasma,  as  origi¬ 
nally  developed  by  Langmuir,  it  is  assumed  that  the  mean-free  path  of  all  plasma  par¬ 
ticles  is  large  in  comparison  with  the  probe  dimensions,  in  order  that  the  plasma 
density,  temperature,  and  potential  as  derived  from  the  measured  characteristics  be 
the  same  as  those  of  the  plasma  undisturbed  by  the  probe.  When  this  condition  does  not 
hold,  the  plasma  density  and  potential  in  the  vicinity  of  the  probe  are  different  from 
those  of  the  undisturbed  plasma.  The  perturbation  of  the  plasma  by  a  probe  is  analyzed 
in  this  report  as  a  problem  in  ambipolar  diffusion,  subject  to  the  assumptions  that  the 
mean-free  path  of  plasma  particles  is:  (a)  comparable  to  or  smaller  than  probe  dimen¬ 
sions,  and  (b)  much  greater  than  the  thickness  of  any  sheaths  that  may  be  present. 
Sheaths  must  therefore  be  assumed  to  be  thin  in  comparison  with  probe  dimensions . 
Analyses  are  carried  through  for  spherical  probes  and  for  plane  probes  oriented  normal 
to  a  magnetic  field. 

The  results  can  be  expressed  in  terms  of  a  parameter  Q  which  is  apprc  lately 
equal  at  zero  sheath  potential  to  the  sum  of  the  ratios  of  probe  size  to  electi  mean- 
free  path  and  probe  size  to  ion  mean-free  path.  Since  Q  also  depends  on  sheath  poten¬ 
tial,  the  current-voltage  characteristics  are  distorted.  Methods  for  determining  zero 
sheath  potential,  and  calculating  the  properties  of  the  undisturbed  plasma  from  the  probe 
curves  are  given. 
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GLOSSARY 


Symbol 


n 


l^e'  '"i 


Te-  Ti 


V  Vi 


Definition 

Density  of  electrons  or  ions  in  neutral  plasma 

Density  of  electrons  or  ions  in  the  unperturbed  plasma 

Density  of  electrons  or  ions  at  the  sheath  boimdary  (i.e.,  at  the  probe) 

Electron,  ion  mobility 

Particle  current  density  of  electrons  or  ions 

Voltage  equivalent  of  electron,  ion  temperatures  in  the  plasma 


Voltage  equivalent  of  electron,  ion  temperatures  at  the  sheath  boundary 
Ambipolar  diffusion  coefficient 
Probe  radius 

Potential  difference  across  the  sheath,  positive  when  probe  is  positive 
with  respect  to  plasma  outside  the  sheath 

Potential  of  the  unperturbed  plasma 

Potential  of  the  plasma  at  the  sheath  boundary  (i.e.,  at  the  probe) 

Vp-V^  perturbation  of  the  plasma  potential  by  the  probe 
2  2 

Equals  1  +  (where  B  is  magnetic  field) 

Equals  1  +  p?B^ 

Equals  1  +  p^p^B^ 

Equals  z/pR 


I.  INTRODUCTION 


In  the  analysis  of  the  current-voltage  characteristics  of  a  probe  in  a  plasma,  as 
originally  developed  by  Langmuir,^  it  is  assumed  that  the  mean-free  path  of  all  plasma 
particles  is  large  in  comparison  with  probe  dimensions.  When  this  is  the  case  the  elec¬ 
tron  density  and  temperature,  and  plasma  potential  derived  from  the  analysis  are  the 
same  as  those  of  the  plasma  in  the  absence  of  the  probe. 

When  this  is  not  the  case  the  plasma  density  and  potential  in  the  vicinity  of  the  probe 
may  be  quite  different  from  those  of  the  unperturbed  plasma,  since  ions  and  electrons 
must  diffuse  from  the  undisturbed  regions  of  the  plasma  to  the  vicinity  of  the  probe .  If 
the  mean-free  path  is  very  short,  the  density  gradients  in  the  plasma  may  be  quite  large, 
and  the  plasma  density  in  the  vicinity  of  the  probe  will  be  substantially  different  from 
that  of  the  tmdisturbed  plasma.  This  situation  is  aggravated  when  a  magnetic  field  is 
present  because  of  the  reduction  of  the  diffusion  coefficient  of  plasma  particles  in  the 
direction  transverse  to  the  magnetic  field. 

The  analysis  that  is  given  herein  was  undertaken  for  the  purpose  of  interpreting  the 
probe  characteristics  in  a  discharge  tube  with  an  axial  magnetic  field.  Similarly  moti- 
vated  studies  have  been  reported  by  Bickerton  and  von  Engel  and  by  Bertotti.  Neither 
of  these  is  applicable  to  the  present  case.  The  first  study  treated  plane  probes  parallel 
to  the  ma.gnetic  field  and  was  restricted  to  relatively  low  magnetic  fields,  while  in  the 
second  it  was  assumed  that  the  longitudinal  mean-free  path  of  all  particles  is  long  in 
comparison  with  any  relevant  dimension.  Neither  of  these  conditions  is  satisfied  in  the 
discharge  in  which  measurements  were  made. 

Because  the  use  of  magnetic  fields  up  to  60,000  gauss,  or  more,  was  contemplated, 
the  configuration  of  a  plane  probe  normal  to  the  magnetic  field  was  chosen.  We  assume 
the  presence  of  a  sheath  between  probe  and  plasma.  Provided  that  the  sheath  is  thin  in 
comparison  with  the  mean-free  path,  the  random  currents  of  ions  and  electrons  across 
the  sheath  from  plasma  to  probe  will  be  unaffected  by  the  presence  of  the  magnetic  field, 
since  they  depend  only  on  the  component  of  velocities  normal  to  the  probe  and  parallel 
to  the  magnetic  field. 

Although  the  probe  current  for  a  given  plasma  density  immediately  outside  the  sheath 
is  not  affected  by  the  magnetic  field,  the  plasma  density  and  potential  themselves  most 
certainly  are.  The  electrons  and  ions  must  diffuse  to  the  sheath  boundaries,  and  their 
diffusion  coefficients  are  strongly  affected  by  the  magnetic  field. 

The  problem  is  one  of  ambipolar  diffusion;  we  must  match  the  diffusion  current  from 
the  plasma  to  the  random  current  across  the  sheath  toward  the  probe.  The  axial  sym¬ 
metry  of  this  problem  dictated  the  use  of  cylindrical  coordinates,  which  has  the  disad¬ 
vantage  of  yielding  answers  in  the  form  of  an  infinite  series  of  Bessel  functions. 

In  setting  up  the  problem,  we  recognized  that  a  similar  physical  situation  exists  in 
the  absence  of  a  magnetic  field  when  mean-free  paths  are  small  in  comparison  with  probe 
dimensions.  Again,  we  have  to  match  an  ambipolar  diffusion  current  from  the  plasma 
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to  the  random  current  across  the  sheath  toward  the  probe.  Since  this  case  can  be  treated 
in  spherical  coordinates  (as  a  one-dimensional  problem),  provides  considerable  insight 
into  the  physical  basis  of  the  problem,  and  yields  answers  in  a  particularly  simple  closed 
form,  it  will  be  treated  in  some  detail  first. 

The  results  of  this  analysis  are  far  from  being  a  trivial  warm-up  exercise,  since 
many  recent  probe  investigations  of  normal  discharges  in  the  absence  of  magnetic  fields 
have  invaded  the  domain  of  probe  sizes  comparable  to  or  greater  than  mean-free  paths. 
The  results  reported  herein  can  be  used  to  estimate  the  degree  of  perturbation  of  the 
plasma  by  the  probes  in  these  experiments,  and  permit  the  correction  of  the  data  for  the 
errors  introduced  thereby. 

We  note  that  the  procedure  of  matching  diffusion  current  to  random  current  is  similar 
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to  the  treatment  by  Allis  and  Buchsbaum  of  the  "pre-sheath."  This  they  classify  as  the 
region  of  the  plasma  (in  which  m  =  n^)  in  which  the  ions  collected  by  the  probe  are  pro¬ 
duced,  and  into  which  there  is  a  penetration  of  electric  field  from  the  probe  to  satisfy 
the  condition  of  continuity  of  ion  current.  Our  work  differs  from  that  of  these  authors 
in  that  they  consider  ion  mean-free  paths  to  be  large  in  comparison  with  all  dimensions, 
and  assume  electron  density  to  be  everywhere  in  equilibrium  with  electrostatic  potential. 

1  consider  the  case  in  which  all  mean-free  paths  are  short  in  comparison  with  the  dimen¬ 
sions  of  the  pre-sheath,  so  that  diffusion  equations  apply.  The  pre-sheath  is  then  the 
region  of  the  plasma  perturbed  by  the  probe.  Moreover,  since  results  are  to  apply  for 
positive  probe  potentials  at  which  electron  currents  may  be  large,  the  assumption  that 
the  electron  density  is  in  thermal  equlibrium  with  electrostatic  potential  cannot  be  made. 

A  calculation  similar  to  that  imdertaken  here  has  been  reported  by  Davydov  and 
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Zmanovskaja.  These  authors  consider  the  case  of  ambipolar  diffusion  to  a  sheath  sur¬ 
rounding  a  spherical  or  cylindrical  probe  in  the  absence  of  a  magnetic  field  for  ion  tem¬ 
peratures  that  are  negligible  in  comparison  with  electron  temperatures.  Boyd^  has 
criticized  their  results  on  the  ground  that  the  boundary  conditions  at  the  sheath  edge  are 
much  more  complex  than  the  abrupt  discontinuity  between  neutral  plasma  and  unipolar 
sheath  employed  by  Davydov  and  Zmanovskaja.  Boyd  distinguishes  between  three 
regions:  (a)  a  xmipolar  sheath  in  which, ions  are  falling  freely  toward  the  probe;  (b)  a 
quasi-neutral  region  in  which  ions  are  being  accelerated  toward  the  probe,  making  a  few 
collisions  with  gas  atoms,  with  velocity  proportional  to  the  square  root  of  electric  field; 
(c)  a  diffusion  region  where  ions  move  by  ambipolar  diffusion  in  a  neutral  plasma. 

In  the  diffusion  region  the  plasma  density  at  the  inner  boundary  is  determined  in  part 
by  the  radius  of  the  inner  boundary;  the  results  of  Davydov  and  Zmanovskaja  agree  with 
those  of  Boyd  when  this  radius  is  not  sensibly  different  from  that  of  the  probe,  that  is, 
when  the  thickness  of  sheaths  (a)  and  (b)  are  small  in  comparison  with  probe  radius. 
When  this  is  not  true,  the  problem  becomes  much  more  complex  because  the  "effective 
probe  radius"  to  be  used  for  the  ambipolar  diffusion  part  of  the  problem  becomes  a  func¬ 
tion  of  the  potential  drop  across  the  sheath,  and  the  more  elaborate  matching  procedure 
of  Boyd  must  be  used. 
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We  note  that  in  high  electron  density  plasmas  the  thicknesses  of  the  regions  (a)  and 
(b)  can  be  small,  especially  in  the  range  of  only  a  few  volts  on  either  side  of  zero  poten¬ 
tial  drop  across  them.  We  shall  therefore  follow  the  procedure  of  Davydov  and 
Zmanovskaja,  including  the  ion  temperature  explicitly,  and  show  that  in  the  domain  of 
validity  of  this  procedure  ^  10  cm  ,  Tp/X  >  0.1)  information  about  the  plasma  ion 

temperature  can  be  deduced  from  the  shape  of  the  probe  characteristic  near  zero  sheath 
potential. 

In  Section  II  we  shall  consider  a  spherical  probe  in  a  spherical  discharge  tube,  fol¬ 
lowing  closely  the  procedure  of  Davydov  and  Zmanovskaja.  In  Section  III,  a  plane  probe 
normal  to  an  axial  magnetic  field  will  be  taken  up.  Section  IV  will  be  devoted  to  a  dis¬ 
cussion  of  the  results,  and  Section  V  to  recently  reported  probe  studies  in  the  light  of 
these  results .  Insofar  as  possible,  the  stress  will  be  placed  on  physics,  and  mathemat¬ 
ical  details  will  be  given  in  the  appendices. 

The  following  assumptions  have  been  made. 

(i)  The  plasmas  are  those  of  diffusion-controlled  discharges,  in  which  there  is  con¬ 
tinuous  ion  production  that  is  proportional  to  electron  density  to  replace  diffusion  losses 
of  electrons  and  ions. 

(ii)  Electrons  have  a  Maxwellian  energy  distribution;  electron  and  ion  temperatures 
are  assumed  to  be  independent  of  position,  and  to  be  unchanged  by  the  presence  of  the 
probe . 

(iii)  The  probe  radius  is  very  much  smaller  than  the  tube  radius  (<1  per  cent).  This 
insures  that  losses  of  electrons  and  ions  to  the  probe  are  small  in  comparison  with  dif¬ 
fusion  losses  to  the  walls;  this  restriction  is  necessary  if  electron  temperature  is  to  be 
unchanged  by  the  presence  of  the  probe . 

(iv)  Ion  mobility  in  the  absence  of  a  magnetic  field  is  negligible  in  comparison  with 
electron  mobility. 

(v)  Probe  radius  is  much  greater  than  sheath  thickness;  or  alternatively,  the  outer 
radius  of  sheath  is  essentially  the  same  as  that  of  the  probe.  Sheath  thickness  is  small 
in  comparison  with  the  mean-free  path. 

(vi)  "Classical"  diffusion  in  a  magnetic  field  is  assumed  to  apply. 
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n.  SPHERICAL  PROBE  IN  A  SPHERICAL  TUBE 


2.1  PERTURBATION  OF  ELECTRON  DENSITY 


We  must  first  set  up  the  ambipolar  diffusion  equation,  which  we  do  without  making 
the  usual  assumption  that  electron  and  ion  currents  are  equal.  We  assume  only  that 
electron  and  ion  densities  are  equal  and  that  div  Fg  =  div  F^^  =  nv^,  the  rate  of  ion  pro¬ 
duction  per  unit  volume.  We  start  with  the  transport  equations 

Fe  =  -(x^VaV-nl 

(1) 

F.  =  np.E  -  fi.V.V-n  J 


Take  the  divergence  of  both  sides,  set  V-  Fg  =  V  -  F^  =  nv.,  and  eliminate  E  between 
the  two. 


nv.  =  -<V  +V.)  V^n 
1  e  1 


nv. 


V  n  +: 


-=  0 


(2) 


Here, 

with  m  «  ^^.  .  This  is  the  same  equation  that  we  would  have  obtained  if  we  had  assumed 
n.  e  _ 

F_  *  F.  as  Schottky  did  in  the  ambipolar  diffusion  theory,  but  we  derive  it  under  far 

6  X 

less  restrictive  assumptions. 

We  must  solve  this  equation  subject  to  boundary  conditions  at  the  probe  of  radius  Tp 
centered  at  r  =  0,  and  at  the  wall  r  =  R.  At  the  probe,  we  have 


F  (r  )  =  -n  .  jV  €  =-niiE-uV 

e'  p'  P  V  ®  ®  drlr^. 

Tit’^p)  =  ""p  V Sir 


eV, 


(3) 


where 


V  /V 
s'  es 


for  V  <  0; 
s 


Ci=l 


for  V  <  0 
s 


-V  /V. 

for  V  >  0;  €.  =  e  ®  for  V„  >  0 

S  I  s 
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The  left-hand  side  of  Eqs.  3  represents  the  current  of  electrons  and  ions  across  the 
sheath  to  the  probe,  as  determined  by  the  potential  difference  Vg  across  the  sheath  and 
the  plasma  density  at  the  sheath  edge,  np.  The  right-hand  side  represents  the  transport 
of  electrons  and  ions  out  of  the  plasma  to  the  sheath  edge  under  the  combined  influence 
of  density  gradients  and  electric  fields. 

Eliminating  E  between  the  two  equations,  we  get 


1  dn 

n  dr 
P 


where 


Q  =  Qg  + 


(4) 


Q.  = 


e  -  Pg(Vg+V.) 


€ 


e 


r 

P 


Pi(V^-tV.) 


e. 

1 


The  solution  of  (2)  in  spherical  coordinates  is 
n=^{A  sin  X  +  B  cos  x}. 


(5) 


where  x  = 


•r.  When  no  probe  is  present,  B  =  0,  since  the  density  cannot  be 


infinite  at  r  =  0;  A  *  n^  and  *  /  ^  must  equal  tt/R  in  order  that  n  =  0  at  the  walls . 

^  a 


When  a  probe  is  present  with  radius  small  enough  that  x^  =  irr^/R  «  1,  the  value  of 

V.  is  not  altered  by  the  presence  of  the  probe  and  x  remains  equal  to  mc/R.  We  deter- 

mine  the  ratio  B/A  from  boundary  condition  (4)  by  equating  — jP L  from  (5)  to  Q/r 

"  I  P  ^ 

and  solving  for  B/A.  It  is  shown  in  Appendix  II  that 


A  ^ 


% 

1  +  Q’ 


(6) 


If  the  probe  radius  is  very  much  smaller  than  the  tube  radius,  the  entire  region  of  the 
plasma  which  is  perturbed  will  be  that  for  which  n  in  the  absence  of  the  probe  is 

SXn  X  011*1  y 

n  =  n  - « - =  n  .  In  substituting  (6)  in  (5)  we  can  regard  cos  x  »  1  and - »  1 

and  obtain  for  the  density  in  the  perturbed  plasma 
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(7) 


The  plasma  density  at  the  sheath  edge,  which  will  be  determined  from  the  Langmuir- 
probe  analysis,  is 


n 

_  _  -Q 

P  ■  1  +  Q  ■ 


(8) 


The  plasma  density  versus  distance  from  the  probe  given  by  Eq.  7  is  shown  in 
Fig.  1.  Note  that  this  is  a  very  slowly  varying  function.  At  a  distance  of  10  probe  radii 
into  the  plasma  the  perturbation  is  still  10  per  cent  as  large  as  it  is  at  the  probe  surface 
We  shall  have  more  to  say  about  this  point  in  Section  IV. 


Fig.  1.  Electron  density  in  the  plasma  as  a  function  of  distance 
from  a  spherical  probe  for  which  Q  =  1 . 


2.2  PERTURBATION  OF  PLASMA  POTENTIAL 


To  calculate  the  electric  field,  we  return  to  the  transport  equations,  take  the  diver¬ 
gence  of  both  sides,  and  solve  for  the  divergence  of  nE,  noting  that  div  *  div  F^ 


Invoking  the  fact  that  and  assuming  that  will  be  comparable  to  or  less 


than  V^,  we  obtain 

V-(nE)  =  -VgV^n. 


(9) 


In  spherical  coordinates, 

1  8  ,Z„  .  „  1  8  /  2  8n\ 

7  8?  “^r)  •  "'^e  7  if) 


A, 

8r 


(r^nE^)  dr  *  -V^ 
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Integrating  and  solving  for  the  electric  field,  we  obtain 


_  _v  1  dn  4.  k:  ^  1 

®r  -  '^elT^+^-Z  n 


(10) 


where 


K  =  nE+V 


dn  1 


e  dr  Jrp' 


and  can  be  evaluated  by  solving  Eq.  3  for  E^,  the  electric  field  at  the  edge  of  the  sheath. 

%  +V  ) 

It  is  shown  in  Appendix  I  that  K  =  — - fjTQ — • 

P  ^ 


„  _  „  1  dn  ^  Qe<Ve^-Vi)  ^  “o 

®r--VeHd?+  -1+Q  ^p;7 


(11) 


In  the  unpertiu'bed  plasma,  n  a  a  constant  n^,  and  E  »  0.  We  can  regard  the  potential 
of  the  unperturbed  plasma  as  a  constant  V^.  We  shall  call  the  potential  of  the  plasma 
at  the  sheath  boimdary  Vp. 


V  =  V 
o  p 


-r 


E  dr 
r 


=  AV 

p  o 


=r 


E,dr. 


Here,  AV  is  the  perturbation  of  plasma  potential  at  the  sheath  boundary  caused  by  the 
presence  of  the  probe. 


AV  V  1  ‘in  1 

=  na?-‘i^+  1+Q  ^p  J 


dr 


‘P  2 


r  1 1- 


x(l+<3), 


The  first  integral  is  equal  to  -V^  In-jj—  =  “V^  In  (1-tQ). 

P 

Q^(V.+V) 

It  is  shown  in  Appendix  I  that  the  second  is  equal  to  —  In  (1+Q).  Therefore 


AV  =  Jin  (1+Q).  (12) 

The  first  term  in  (12)  is  the  Boltzmann  term,  resulting  from  the  fact  that  the  potential 
and  electron  density  are  related  in  thermal  equilibrium  according  to  n  ~  exp  (V/V^). 
From  this  term  the  plasma  potential  is  shifted  negative  by  the  presence  of  the  probe. 
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The  second  term  is  large  only  when  is  a  large  fraction  of  Q.  Reference  to  Eq.  4 
shows  that  this  can  happen  only  when  the  sheath  potential  is  slightly  negative,  zero  or 
positive.  Then  the  electron  current  to  the  probe  is  large  and  the  electric  field  in  the 
plasma  must  be  reduced  from  the  thermal  equilibrium  value  to  permit  the  increased 
electron  current  to  flow.  Whether  or  not  the  electric  field  actually  reverses,  and  the 
perturbation  of  plasma  potential  becomes  positive,  depends  on  the  exact  value  of  the 
ratios  of  Pj/pg  and  Vj/Vg. 
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m.  DISC  PROBE  NORMAL  TO  A  MAGNETIC  FIELD 


3.1  PERTURBATION  OF  PLASMA  DENSITY  BY  THE  PROBE 


The  basic  procedure  is  the  same  as  before;  set  up  the  ambipolar  diffusion  equation 
under  the  assumption  that  I\  ^  div  F-  =  div  F^  =  nvj^.  Solve  it  subject  to  the  bound¬ 
ary  condition  that  diffusion  current  toward  the  sheath  boundary  from  the  plasma  is  equal 
to  the  random  current  across  the  sheath  to  the  probe.  This  determines  electron  density 
as  a  function  of  position;  then  use  the  transport  equations  to  calculate  the  electric  field 
that  is  necessary  to  bring  about  the  flow  of  a  given  probe  current. 

The  derivation  of  the  ambipolar  diffusion  equation  in  the  presence  of  a  magnetic  field 
is  somewhat  more  complicated.  Allis  and  Buchsbaum  show  that  when  the  azimuthal 
components  of  current  flow  are  eliminated  the  transport  equations  in  cylindrical  geom¬ 
etry  are: 


dn  N 


Etr  vAii 

—  u.  V  -3 — 

'  e  z  f^e  e  dz 

nuL  E  u.  V  . 
^e  r  ^e  e  dn 


er 


2 

Ye 


2  dr 
Ye 

dn 


F.  =  nu..E  -  u.V.  ^ 
iz  '^1  z  1  dz 


^  ir  2  2  dr 

Yi  Vi 


(13) 


Here,  =  1  +  y?  =  1  +  p?B^.  With  only  the  assumption  that  div  F^  =  div  F^  =  nv^ 

there  are  too  few  equations  to  eliminate  all  components  of  E  and  F.  I  choose  to  elim¬ 
inate,  for  the  present,  both  components  of  the  electric  field.  This  procedure  gives  the 
result 


—  F  +  — F.  =-(V+V.)-^ 
u  ez  u.  iz  '  e  1  dz 


(14) 


Yi 


dn 


—  +-^r,^  =  -(v+v,)^ 


er  n- 


ir 


e”i'dr  J 


Note  that 


1  a  dF  dF 

—  |-(rF  )  =  div  F  '  =  nv.  -  — 5-^ 

r  8r  er  e  dz  x  dz 


- 1-  (rF.  )  =  div  F.  - 
r  3r  '  ir  1 


dF. 

iz 


.  -  nv.  . 

dz  1  dz 


(15) 


By  taking  the  z  derivative  of  the  z  component  of  Eq.  14  and  operating  with  the 
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r  component  of  the  divergence  on  the  r  equation,  with  the  aid  of  (15),  we  can  eliminate 
all  but  one  of  the  components  of  the  divergence  of  r,  or  r„.  I  choose  to  retain  dr  /dz. 

1  c  0Z' 

straightforward  differentiation  and  algebraic  manipulation  leads  to 


The  left-hand  side  of  ( 16)  is  a  normal  homogeneous  diffusion  equation  with  all  of  the 
usual  solutions  obtained  by  setting  it  equal  to  zero.  Because  of  the  right-hand  side, 
there  is  also  a  particular  integral,  a  function  of  z  and  r,  which  leads  to  a  whole  new 
family  of  solutions.  We  restrict  ourselves  to  four  limiting  cases: 

dr 

(i)  -■  =  0.  The  right-hand  side  of  Eq.  16  is  zero  and  it  becomes 


p2  r  8r  ar  ) 


Y- 

+  ^+nv. 


0, 


(17) 


2  2 

where  D  =  a.(V  +V.),  p  =  1  +  p.  p.B  ,  and  we  have  assumed  u,  »  p,.  throughout.  This 
case  would  represent  conditions  near  the  anode  of  a  discharge  tube  carrying  a  current 

that  is  independent  of  z  in  an  axial  magnetic  field. 

dr.  dr  p  .2 

aZ^  =  From  (14)  we  have  ^  -  jig(V^+Vj^)'2-§. which  wdien  substituted 


dz 

in  (16)  gives 


dz 


i.  (r^\ 
p2  r  ar  V  dr/ 


+  D  ^  ^  +  nvi  =  0. 
^  dz^  ^ 


(18) 


This  equation  was  previously  derived  by  Allis'*  under  the  assumption  that  E  =0. 

dr  dr 

(iii)  — =  nv^.  Since  div  r^  '*“r  ^  ^^^er^  “  ”’'i*  assumption  is  tanta¬ 

mount  to  assuming  that  the  radial  component  of  the  divergence  of  the  electron  current 
flow  is  negligible  in  comparison  with  both  the  z  component  and  the  production  per  unit 
volume. 


Under  this  assumption,  Eq.  16  becomes 


2 

Yi 


r  ar  V  drj 


+  D. 


a^n 


+  nv.  =  0. 


(19) 


•7 

Equation  19  was  previously  derived  by  Simon. 

^^ez  1  a 

(iv)  =  “  7"  aP  "’'i'  terms  on  the  right-hand  side  being 

of  the  same  order  of  magnitude.  It  is  shown  in  Appendix  II  that  this  assumption 
leads  to 


1  8 


.2  r  ar 
P 


+  D 


(20) 
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This  is  the  "normal"  ambipolar  diffusion  equation  that  would  be  obtained  by  letting 

^i  =  r,- 

The  question  now  arises.  Which,  if  any,  of  these  four  diffusion  equations  is  appli¬ 
cable  to  the  present  case?  We  can  rule  out  the  first  two  on  the  ground  that  neither 

dr  dr. 

©z  xz  ■ 

— nor  -  can  be  zero  for  the  plane  probe  normal  to  the  z-direction.  We  are  left 

with  the  choice  between  the  Simon  and  the  normal  ambipolar  equations .  Here,  the  choice 

1  8 

rests  squarely  on  the  relative  magnitudes  of  nv^^  and  —  gp  (rr  ). 

For  the  infinitely  long  positive  column  in  a  magnetic  field  with  diffusion  to  the  walls 
the  only  loss  mechanism,  these  two  are  equal.  The  Simon  diffusion  equation  was  origi¬ 
nally  derived  for  a  relatively  short  discharge  column  in  an  axial  magnetic  field  with 
electron  loss  to  conducting  end  plates  being  the  major  loss  mechanism.  For  this  case 
1  8 

—  ^  negligible  in  comparison  with  and  nv^,  and  it  becomes  more  and 

more  negligible  the  higher  the  magnetic  field. 

We  have  from  the  outset  restricted  ourselves  to  considering  probes  that  are  small 
in  radius  in  comparison  with  tube  radii.  If  the  probe  radius  is  sufficiently  small  that 
electron  loss  to  the  probe,  even  though  it  may  occur  from  a  long  column  of  plasma  of 
radius  approximately  equal  to  probe  radius,  shall  be  small  in  comparison  with  the  radial 
diffusion  loss  to  the  tube  walls  over  the  "characteristic  length"  of  the  positive  column, 
then  the  electron  temperature  and  ionization  frequency  in  the  plasma  will  be  unchanged 
by  the  presence  of  the  probe.  We  therefore  employ  the  normal  ambipolar  diffusion 
equation. 

It  is  convenient  to  express  this  in  dimensionless  variables,  p  =  -y,  t  =  -w-,  and  let 

(D  n  pK 

~  •  Equation  20  then  becomes 
2 

18/  8n\  ■  8  n  .  ,.2  „ 

Equation  21  is  to  be  solved  subject  to  the  boundary  conditions: 

for  4  =  00  n  =  n^J^(2.4p) 

At  i  =  0,  at  the  probe,  for  r  <  rp, 

/5I  .  ,  .  vin  'l 

~”p  V  2irm  ^e  ~  "^p^e  z  ^^e  e  dz  „ 
e  z  =  0 


,  e  -  n  IT  _  V  ^ 

^p  Y  2irm.  ^i  "  "p^'^i^z  ^‘i^i  dz 

X  Z  s 


Equations  22  are  essentially  the  same  as  Eq.  3  and  lead  to  the  same  bound¬ 
ary  condition: 
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(23) 


1  dn 

ir-ai 

p 


Q 

r 


for  r  <  r  , 
P 


where  Q  is  as  we  have  defined  it  in  Eq.  4.  Since  the  whole  problem  is  symmetrical 

dn 

about  4  =  0,  in  the  plasma  for  r  >  rp  we  must  have  =  0  at  C  =  0*  otherwise  the  density 
gradient  would  have  a  discontinuity  at  4  =  0.  Substituting  p  for  r  in  (23),  we  can  sum¬ 
marize  the  boundary  conditions: 

dn 


t  =  0 


P  <  Pr 


P  >  Pr 


(24) 


i  =  00 

It  is  shown  in  Appendix  II  that  the  solution  of  (21),  subject  to  conditions  (24),  is 


00 

n(p,4)  =  n^J^(2.4p)  -2pQnp  ^ 


exp  (-nAJ-(2-4)^c)  Ji(m^Pp)  J„(m.p) 


i=2 


m 


3 


^m^-(2.4)^  J^(m,) 


J 


•th 


Here,  m^  is  the  j  zero  of  J^,  the  Bessel  function  of  zero  order.  To  determine  the 
ratio  of  np  to  n^,  set  p  =  0,  4  =  0-  The  left-hand  side  then  becomes  np,  the  plasma 
density  at  the  sheath  edge,  by  definition.  Solve  for  n 


1  -f  2pQ 


(25) 


j=2  m. 


The  series  is  a  series  of  alternating  terms  of  decreasing  magnitude,  and  therefore  con¬ 
verges,  but  the  convergence  is  very  slow.  It  can  be  greatly  improved  by  calculating 
n(p=0,  ^=e),  where  €  is  small  in  comparison  with  the  sheath  thickness.  This  will  differ 
by  a  negligible  amount  from  the  plasma  density  at  the  sheath  boundary;  the  retention  of 
the  decreasing  exponential  terms  in  the  series  greatly  improves  the  convergence,  and 
in  fact,  makes  it  absolute. 

Note,  however,  that  in  the  limit  of  zero  magnetic  field,  p  =  1,  and  we  are  dealing 
with  the  same  problem  that  we  have  treated  in  Section  II,  although  in  a  different  geom¬ 
etry.  Equation  7  shows  that  the  perturbation  of  plasma  density  extends  many  probe  radii 
into  the  plasma.  Figure  2  is  intended  to  suggest  that  outside  of  a  certain  radius,  which 
is  somewhat  greater  than  the  probe  radius,  the  contours  of  constant  plasma  density  will 
be  much  the  same  whether  the  probe  is  spherical  or  plane. 

Therefore,  physical  intuition  tells  us  that  the  perturbation  of  the  plasma  by  the  probe 
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Fig.  2.  A:  contours  of  constant  electron  density  in  the  plasma  as  a  function  of  dis¬ 
tance  from  a  plane  probe.  B:  contours  of  constant  electron  density  in  the 
plasma  as  a  function  of  distance  from  a  spherical  probe.  The  difference 
between  A  and  B  is  slight  except  very  near  the  probe. 


will  be  essentially  the  same  whether  it  is  a  plane  probe  or  a  spherical  one  of  the  same 
radius,  within  numerical  factors  of  the  order  of  unity.  We  conclude,  then,  that  the 
infinite  series  in  Eq.  25,  which  does  not  depend  on  magnetic  field,  has  a  sum  of  approx¬ 
imately  1/2.  For  the  plane  probe  in  a  magnetic  field,  therefore,  Eq.  8  is  replaced  by 


“  1  +  PQ 


(26) 


The  effect  of  the  magnetic  field  then  is  to  magnify  the  Q-value  for  a  particular  probe 

radius  —  mean-free  path  combination  by  the  factor  p  =  \J  1  +  p  p-p^-  In  order-of- 
magnitude  numbers,  this  is  a  factor  of  10  at  1  mm  pressure  and  10,000  gauss. 

3.2  PERTURBATION  OF  PLASMA  POTENTIAL  BY  THE  PROBE 


The  procedure  is  the  same  as  that  in  Section  II.  Calculate  the  electric  field  from 
the  transport  equations,  and  integrate  from  the  probe  to  infinity  to  find  the  difference 
in  potential  between  the  plasma  at  the  probe  sheath  surface  and  the  undisturbed  plasma 
at  infinity. 

We  start  with  Eqs.  13,  and  obtain  the  divergence  of  and  F^,  setting  both  equal  to  nv^ 


<Uv  r,  .  nv,  .  ±  („E^)  -  ^  — I  1  i  (rnEp  — ^  7  ^  (--If) 

q  •  "“1  ■  A  7  h  -^7  S  • 
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1  8 


Eliminate  —  ^  (rnE  )  between  the  two  equations, 
r  or  r 


/  2 

.  Ve  / 

2  2\ 

Ci 

< 

il 

Vi-^ej 

(27) 


2  2  ^  2„2  ^  .  2  _2  ,  „2\ 
\  ^e  ..  ^  ^  ^  ►^i  ^  >^eV^>^e>^i^  ) 


t*i 


=  <VVe)^- 

a 

Here,  we  have  neglected  in  comparison  with  |j.^.  From  Eq.  20,  we  have 

(■•If)  - 

Substituting  in  Eq.  27  and  solving  for  E^  gives 


8z  <^z) 


Ye-n 


8z 


Straightforward  algebraic  manipulation  shows  that  this  reduces  to 


a  .  „  ,  d^n 


(28) 


This  is  basically  the  same  equation  as  the  one  that  we  obtained  for  the  spherical  case, 
except  that  here  we  have  z  derivatives  instead  of  r  derivatives.  The  fact  that  the  mag¬ 
netic  field  drops  out  as  far  as  the  relationship  between  E^  and  dn/dz  is  concerned  is 
one  that  we  could  have  deduced  by  intuition,  since  plasma  particles  diffuse  freely  par¬ 
allel  to  B. 

Neglecting  in  comparison  with  and  integrating  (28)  gives 


E  =-V  + 

z  e  n  8z  n 

We  determine  K(r)  from  the  value  of  E  ,  n  ,  and  dn/dz  at  rsO,  z^O. 

z  p 


dn 


K(0)  =  n  E  +  V 
'  '  p  p  e  dz 


(29) 


Since  the  boundary  conditions  for  the  plane  probe  at  zsO  are  the  same  as  those  for  the 
spherical  probe  at  rsr^,  K  has  the  same  value  as  that  derived  in  Appendix  I: 
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npQe(V,+V.) 


K(0)  =- 


In  this  case,  however,  n  =  .  ,  -  ,  so  that  (29)  becomes 

P  i  +  pw 

,  ,  n  Q  (V  +V.) 
p, ^  1  dn  o  e  e  i' 

~  *e  n  dz  (l+pQ)rp-n 


(30) 


where  this  applies  only  for  r=0. 

Note  that  the  exact  integral  of  this  expression  from  z=0  to  z=m  along  r=0  will  yield 
the  perturbation  of  plasma  potential.  This  will  involve  integrating  over  infinite  series, 
which  we  do  not  consider  worth  while  now.  We  shall,  therefore,  resort  to  an  approx¬ 
imation  for  the  variation  of  n  with  z,  as  shown  in  Fig.  3. 

With  this  approximation 


n 


1  +  PQ 


=  n 


z  3^  pr 


and  we  carry  out  the  integration  between  the  limits  of  np  and  n^  (z=0  to  z=prp): 


-  Z/Pfp 


Fig.  3. 


A:  Schematic  representation  of  correct  variation  of  electron  density  vs 
distance  along  the  z  axis  away  from  a  plane  probe.  B;  The  function  n^/n^s 

pr^^(l+Q  /r  )  used  for  the  integration  of  Eq.  30  in  Section  III.  C:  The 
1  T  Z  P 


fimction  *ll 

in  Appendix  II. 


PQ 

TTpQ 


exp(-z/prp) 


used  for  the  integration  of  Eq.  30 
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1 


1  dn  .  Qe<Ve+Vi) 


~  ~^e  n  3z  ^  r  1  +  Qz/r 


AV  =  -V„  In 
e 


r"p  _ 

Jo  1  + 


Qz/r 


Q-(V  +V  ) 

AV  =  -V„  In  (1+pQ)  +  — -  In  (1+pQ) 

6  Vi 


AV=-  V^- 


Qe(V^+V.) 


Q 


In  (1+PQ) 


(31) 


In  the  limit  of  zero  magnetic  field  this  is  seen  to  reduce  to  the  same  value  as  for  the 
spherical  case. 

It  may  be  objected  that  the  electron  density  function  chosen  in  Fig.  3  is  highly  arti¬ 
ficial,  but  it  is  shown  in  Appendix  III  that  the  function 


n  =  n  1  - 


PQ  -z/pr 


1  +  PQ 


(dotted  line  in  Fig.  3)  for  the  second  term  in  the  integral  yields 

PQe(V-+Vi) 

in  (1+PQ) 


which  is  sensibly  different  from  Eq.  3 1  only  for  Q  «  1,  for  which  a  diffusion  theory  will 
not  apply. 
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IV.  DISCUSSION 


From  a  perusal  of  the  equations  derived  in  Sections  II  and  III  it  is  clear  that  the  ion 
temperature,  both  in  the  plasma  and  at  the  sheath  boundary,  affects  the  results  through 
its  appearance  in  the  Q's.  In  low-pressure  plasmas  in  which  Vg»V.,  the  ion  temperature 
does  not  affect  Q  .  However,  the  presence  of  V.  (the  ion  temperature  at  the  sheath 

G  IS 

boundary)  in  the  exponential  of  Q.  shows  that  this  quantity  will  have  a  major  effect  upon 
the  exact  shape  of  the  probe  characteristic  for  positive  sheath  voltages . 

The  question  therefore  arises:  What  is  the  value  of  ion  temperature  at  the  sheath 
boundary?  Since  the  ions  are  accelerated  toward  the  sheath  by  the  ambipolar  electric 
fields,  they  might  be  expected  to  arrive  at  the  sheath  edge  with  a  higher  energy  than 
they  possess  in  the  plasma,  Tonks,  in  fact,  derived  the  relationship,  for  long  mean- 
free  paths  and  no  collisions,  that  the  ions  arrive  at  the  sheath  boundary  with  an 
energy  approximately  equal  to  the  electron  mean  kinetic  energy.  Allis  and  Buchsbaum^ 
also  derive  this  result  in  their  treatment  of  the  pre -sheath. 

These  workers  have  all  been  interested  in  the  ion  saturation  current  for  very  nega¬ 
tive  probe  potentials.  The  only  effect  of  the  ion  temperature  at  the  sheath  boundary  is 
to  determine  the  ion  current  density.  The  effect  of  ion  energy  on  the  probe  character¬ 
istic  at  positive  probe  potentials  has  been  ignored.  There  have  been  sound  reasons  for 
this,  of  course;  after  all,  the  electron  current  at  zero  sheath  potential  is  several  hun¬ 
dred  times  the  ion  current,  and  the  decrease  of  ion  current  to  zero  with  positive  sheath 
potentials  cannot  even  be  noticed. 

In  the  regime  to  which  the  results  of  this  study  apply,  there  is  a  much  greater  effect 
of  ion  current  on  the  characteristic  for  positive  sheath  potentials.  As  sheath  voltage 
becomes  positive  and  ions  are  repelled  from  the  probe,  -*0.  Since  is  of  the  same 
order  of  magnitude  as  Q  ,  the  value  of  Q  decreases  substantially  as  Q.  ••0.  Therefore 
the  plasma  density  at  the  sheath  boundary,  n^,  increases  considerable  as  Q^-*0.  This 
is  necessarily  reflected  as  an  increase  in  electron  current  to  the  probe  as  the  ions  are 
repelled  from  the  probe. 

Therefore,  although  the  change  in  ion  current  is  much  too  small  to  measure  directly, 
the  change  in  electron  current,  because  of  the  retarding  of  the  ions,  may  be  quite  large  — 
a  factor  of  2,  or  more.  Moreover,  as  will  be  made  evident  presently,  the  shape  of  the 
characteristic  for  positive  sheath  potentials  is  considerably  affected  by  the  ion  temper¬ 
ature  at  the  sheath  boundary. 

We  are  unable  to  use  the  result  of  Jonks;  the  distance  over  which  the  ambipolar 
electric  field  extends,  and  over  which  the  ions  are  accelerated  toward  the  probe,  is 
several  probe  radii;  that  is,  large  in  comparison  with  the  ion  mean-free  path.  The  ions 
are  therefore  making  many  collisions  on  their  way  to  the  probe  and  losing  energy  in 
large  chunks. 

In  fact,  the  energy  with  which  the  ions  arrive  at  the  sheath  boundary  is  just  equal 
to  the  energy  that  they  had  after  their  last  collision  plus  the  energy  gained  from  the 
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accelerating  field  in  the  last  free  flight,  or 


where  is  the  ambipolar  electric  field  at  the  sheath  edge,  is  the  gas  temperature, 
and  is  the  ion  mean-free  path.  The  factor  of  2  holds  only  for  equal  mass  of  ion  and 
gas  atom.  It  is  shown  in  Appendix  III  that  this  leads  to  the  conclusion  that 

Figures  4,  5,  6,  and  7  show  calculated  probe  curves  for  various  gas  pressure  times 
probe  radius  products,  calculated  for  mercury  ions  in  argon,  under  the  assumption  that 
V.  =  1/30,  l/lO,  l/3,  1  multiplied  by  V  ,  respectively.  These  curves  apply  only  to 
probes  for  which  sheath  thickness  is  small  in  comparison  with  probe  dimension,  for 
which  a  sharp  "break"  at  plasma  potential  is  expected  in  the  unperturbed  case. 

Note  that  all  of  the  curves  approach  the  straight  line  at  negative  probe  potentials; 
thus  it  is  indicated  that  the  measurement  of  electron  temperature  can  still  be  made .  The 
deviation  of  the  curves  from  the  Maxwellian  line  as  zero  sheath  potential  is  approached 
is  the  result  of  the  second  term  in  Eq.  12.  This  is  the  shift  in  potential  in  the  positive 
direction  (from  the  thermal  equilibrium  potential)  required  to  draw  electrons  in  toward 
the  probe. 

The  most  surprising  feature  of  these  curves  is  the  fact  that  the  "knee"  occurs  not  at 
zero  sheath  potential  but  at  a  positive  sheath  potential.  The  reason  for  this  is  that  elec¬ 
tron  current  continues  to  increase  with  positive  sheath  potentials  as  ions  are  repelled, 
Q^-*0  and  the  electron  density  at  the  sheath  edge  increases  from  l/l+Q^+Q^^  to  l/l+Q^. 
Saturation  electron  current  is  not  achieved  essentially  until  all  of  the  ions  are  repelled 
and  electron  density  at  the  sheath  boimdary  no  longer  changes  with  increasing  positive 
sheath  potential. 

We  are  led  to  the  conclusion  that  there  is  no  clear-cut  experimental  method  of  identi¬ 
fying  zero  sheath  potential  without  some  additional  information.  At  high  ion  temper¬ 
atures,  the  "knee"  disappears  completely,  while  at  lower  ion  temperatures  the  positive 
displacement  of  the  "knee"  from  zero  sheath  potential  depends  on  the  ion  temperature. 

It  is  necessary  to  resort  to  a  series  of  successive  approximations  and  calculated  values 
of  the  Q's  to  be  able  to  Identify  zero  sheath  potential  from  a  given  set  of  experimental 
data. 

The  Q's  cannot  be  calculated  without  knowledge  of  electron  and  ion  temperatures. 
The  electron  temperature  can  be  obtained  in  the  usual  way,  from  the  slope  of  the  straight 
line  portion  of  the  logarithmic  plot  of  electron  current  against  probe  potential.  In  order 
to  estimate  the  ion  temperature,  the  following  empirical  procedure  can  be  used.  Meas¬ 
ure  the  potential  difference  between  the  extrapolated  electron-current  vs  probe -potential 
straight  line,  and  the  "knee"  of  the  experimental  curve.  Since  the  "knee"  may  be  diffi¬ 
cult  to  identify  and  the  procedure  is  empirical,  it  is  preferable  to  measure  this  potential 
difference  at  some  lower  current  than  saturation,  say  90  per  cent  of  saturation  (see 
Fig.  8).  Figure  9  shows  this  potential  difference,  AVj^  vs  V^/V^  for  various  values  of 
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Fig.  4.  Probe  electron  current  vs  probe-to-plasma  potential  difference  for  =  30  V. 

and  various  values  of  Q  (7^=0).  Arrows  indicate  location  of  zero  sheath 

potential.  Pressure  and  probe  radius  products  giving  those  values  of  Q  for 
mercury  ions  in  argon  with  1  ev  electron  temperature  are  also  given. 
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Fig.  5.  Probe  electron  current  vs  probe -to -plasma  potential  for  =  10  V^,  and  the 

same  values  of  pressure  times  probe  radius  as  in  Fig.  4.  Arrows  indicate 
'  location  of  zero  sheath  potential. 
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V  =3V. 


Q 


Fig.  6.  Probe  electron  current  vs  probe -to-plasma  potential  for 
Vg  =  3V^,  and  the  same  values  of  pressure  times  probe 

radius  as  in  Fig.  4.  Arrows  indicate  location  of  zero 
sheath  potential. 
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same  values  of  pressure  times  probe  radius  as  in  Fig.  4.  Arrows  indicate 
location  of  zero  sheath  potential. 


pressure  times  probe  radius  p  •  r^,  under  the  assumption  that  there  are  mercury  ions 
in  argon. 

From  a  comparison  of  an  experimental  value  of  with  calculated  curves  of  vs 
V^/Vg  (obtained  from  calculated  curves  of  electron  current  vs  probe  potential)  the  ratio 
V./V  can  be  determined.  Hence  the  value  of  Q  and  Q.  can  be  calculated. 

1  ©  6  X 

Next,  using  the  calculated  values  of  Q  (V  =0),  Q.  (V  =0),  we  determine  the  electron 

6  S  1  o 

current  for  which  the  potential  difference  between  the  extrapolated  i^  vs  line  and  the 
experimental  curve  is 


Fig.  8.  Construction  for  obtaining  A,  measured  electron-current  vs 

probe-voltage  characteristic.  B,  extrapolated  linear  portion  of  the 
characteristic.  Draw  a  line  at  a  current  90  per  cent  of  the  satu¬ 
ration  current;  is  the  potential  difference  between  the  inter¬ 

sections  of  this  line  with  A  and  B. 


23 


5 


pr  =  0.720  MM-CM 
P 


i 


Fig.  9.  AVjj,/Vg  vs  Vg/Vj  for  various  values  of  T  Tp  for  mercury  ions  in  argon. 


(See  Fig.  10.)  This  is  the  electron  current  for  Vg=0.  Then,  extrapolate  the  i^  vs  Vp 
line  for  a  potential  difference,  AV2.  beyond  this  current. 


AV2  = 


V  in  (1+Q^+Q.) 


V^  =  0 


The  current  at  this  potential  is  the  random  current  density  in  the  undisturbed  plasma. 
The  probe  potential  at  this  point  is  the  potential  of  the  undisturbed  plasma,  and  the 
potential  difference  between  this  point  and  the  point  of  zero  sheath  potential  is  given  by 
Eq.  12; 


AV 


In  (1+Qg+Qj) 
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PROBE  POTENTIAL 


Fig.  10. 


Construction  for  determining  zero  sheath  potential.  A  and  B  as  in  Fig.  8. 

Q  (V  +V.) 

Find  the  current  for  which  AV,  =  ^  p.  ^  In  (1-fQ  +Q.)  .  Intersection 

1  '<e  ®  ^  V  =0 

of  this  current  with  A  is  at  zero  sheath  potential  Unperturbed  plasma 
potential  is  AV2  volts  positive  from  the  intersection  of  this  current  with  B, 

where  AV-  =  V  In  (1+Q  +Q.)|  .  Current  at  which  B  intersects  unper- 

^  e  e  1  y  _Q 


turbed  plasma  potential  is  the  random  current  to  a  probe  of  the  same  area 
in  an  unperturbed  plasma.  As  a  check  on  consistency,  the  currents  i  .  i  , 

11  os 

i  should  be  in  the  ratio  i  :  i„:  i  =  1;  .  .  p.  :  .  ■  r»  j.  PT 
p  o  s  p  1  +  Q-  1  +  Q.  +  Q. 
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As  a  check  on  the  correctness  of  the  results,  the  three  currents  i  ,  i  ,  i  ought  to  be 

o  s  p 

in  the  ratio 

....  1  1 

^o=^s=v “  =1 +Q  +Q. 

^  ®  ^  V  =  0 

s 

If  they  are  not,  the  most  likely  source  of  error  is  in  the  determination  of  i  ,  and  there- 

s 

fore  of  0.9 i_,  and  thus  of  V..  Especially  at  high  ion  temperatures,  the  curve  will  not 
reach  saturation  until  relatively  high  positive  sheath  potentials  are  reached,  and  ioni¬ 
zation  by  electrons  in  the  sheath  will  cause  the  measured  current  to  increase  with 
increasing  probe  potential,  and  the  experimental  curve  never  flattens  off.  Therefore, 

make  an  appropriate  adjustment  in  the  assumed  value  of  i  ,  redetermine  the  ion  temper- 

s 

ature,  recalculate  the  Q's,  relocate  zero  sheath  potential,  and  repeat  the  entire  process 
until  it  is  self-consistent. 
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V.  ANALYSIS  OF  PREVIOUSLY  PUBLISHED  WORK 


The  purpose  of  this  section  is  to  examine  a  limited  number  of  recently  published 
high-caliber  probe  studies  in  the  light  of  the  equations  derived  in  this  report.  No 
attempt  is  made  to  make  a  complete  survey,  nor  is  there  any  intention  to  deprecate  fine 
pieces  of  work. 

The  works  considered  are  those  of  Anderson,^  Bills,  Holt,  and  McClure,^®  Medicus,^^ 
12  13 

Verweij,  and  Waymouth.  Of  these,  Medicus  employed  a  shperical  probe  and 
Waymouth  a  plane  probe.  The  others  used  wire  probes.  Since  wire  probes  have  not 
been  analyzed  in  the  present  study,  we  shall  first  consider  the  definition  of  an  "effective 
radius"  for  a  wire  probe.  This  may  be  defined  as  the  radius  of  the  spherical  probe  that 
causes  the  same  perturbation  in  plasma  density. 

It  seems  reasonable,  in  the  light  of  the  physical  mechanism  that  has  been  considered, 
that  two  probes  that  draw  the  same  current  from  a  given  plasma  will  cause  the  same 
degree  of  perturbation.  Since  probes  draw  current  in  approximate  proportion  to  their 
areas,  we  can  say  that  the  cylindrical  probe  has  an  effective  radius  equal  to  the  radius 
of  the  sphere  with  the  same  surface  area: 


r 


eff  " 


where  rp  is  the  radius,  and  I  is  the  length  of  the  cylindrical  wire  probe.  Undoubtedly, 
this  is  only  approximately  true,  but  to  improve  on  it  would  require  a  good  deal  more 
calculation  than  seems  to  be  warranted  now. 

A  second  point  to  consider  in  discussing  the  reported  probe  studies  is  the  fact  that 
probe  lead-in  wires,  or  supports,  are  always  shielded  to  prevent  collection  of  electrons 
or  ions.  Since  shields  are  frequently  of  substantially  larger  radius  than  the  probe  itself, 
especially  for  wire  probes,  diffusion  to  the  shield  may  reduce  the  plasma  density  sub¬ 
stantially  in  the  vicinity  of  the  shield.  This  is  particularly  a  problem  with  wire  probes, 
for  which  a  substantial  fraction  of  the  probe  area  may  lie  within  the  perturbed  region. 

Customarily,  the  probe  shield  operates  at  floating  potential;  the  value  of  Q  for  it  is 
therefore  equal  to  Q^.  By  calculating  Q.  for  the  end  of  the  shield  and  regarding  it  as  a 
spherical  probe,  we  can  use  Eq.  7  to  find  the  variation  of  plasma  density  as  a  function 
of  distance  from  the  shield  (see  Fig.  11).  We  can  therefore  integrate  Eq.  7  to  find  the 
average  electron  density  seen  by  the  probe,  assuming  that  the  Q  for  the  probe  itself  is 
small.  For  a  wire  probe,  this  is  given  by 


Q. 


where  £  is  the  probe  length,  r^j^  the  shield  radius,  and  the  of  the  shield.  It 
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Not  given  by  author. 


GLASS  PROBE  SHIELD 


CYLINDRICAL  WIRE  PROBE 


z 


Fig.  11.  Electron  density  in  the  plasma  vs  distance  from  the  probe  shield  for  a 
cylindrical  wire  probe.  Q  for  the  probe  itself  is  assumed  negligible, 
and  Q.  is  unity. 


is  plain  that  for  probes  that  are  long  in  comparison  with  shield  radius,  the  value  of  n^ 
approaches  n^. 

When  the  Q's  of  probe  and  shield  are  comparable,  and  the  probe  is  not  long  in  com¬ 
parison  with  the  shield  radius,  the  whole  diffusion  problem  becomes  much  more  compli¬ 
cated  and  the  calculation  of  the  corrections  to  be  used  with  wire  probes  is  dubious. 
Finally,  for  the  majority  of  these  cases,  assumption  5  (Section  I)  is  not  really  satisfied, 
and  the  corrections  calculated  here  can  be  regarded  as  valid  to  order  of  magnitude  only. 

Table  I  shows  various  pertinent  data  from  the  authors  cited,  together  with  the  calcu¬ 
lated  values  of  the  Q's. 

It  is  plain  that  of  the  five  references  cited,  only  the  work  of  Medicus  and  of  Bills, 
Holt,  and  McClure  appears  to  avoid  density  and  potential  perturbations.  It  must  be 
pointed  out  that  a  diffusion  theory  is  not  expected  to  apply  when  either  of  the  Q's  is  less 
than  0.1.  Since  Q^,  Q^  are  approximately  equal  to  the  ratio  of  electron  or  ion  mean- 
free  path  to  probe  radius,  and  the  perturbation  in  density  extends  approximately  10  probe 
radii  into  the  plasma,  for  Q^,  Q^  <  0.1,  the  entire  perturbation  would  take  place  in  a 
distance  of  one  mean-free  path  or  less. 

Varweij  has  employed  probe  measurements  to  measure  electron  density,  axial  elec¬ 
tric  field,  and  electron  temperature  in  a  discharge  in  mercury  vapor  and  argon.  Figure  12 
shows  his  published  data  for  electron  density  as  a  function  of  argon  pressure,  together 
with  a  corrected  value  that  I  have  estimated.  Since  Verweij's  measurements  apply  to 
the  case  for  which  «  V  ,  a  well-defined  saturation  current  would  be  expected,  and  the 
saturation  current  would  yield  =  n^/(l+Q^).  It  is  also  possible  to  estimate  too  high 
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(MM) 

Fig.  12.  Data  from  Verweij^  ^  on  electron  density  and  temperature  in  a  discharge  at 
400  ma  in  mercury  vapor  at  7  microns  plus  argon  vs  argon  pressure  (experi> 
mental  points),  together  with  estimated  corrections  for  plasma  perturbation 
(solid  lines). 

an  electron  temperature,  if  the  "best  possible  straight  line"  is  drawn  through  the  experi¬ 
mental  data  right  up  to  the  apparent  zero  of  sheath  potential,  instead  of  relying  only  on 
the  very  negative  sheath  voltage  range.  Figure  12  also  includes  Verweij's  data  on  elec¬ 
tron  temperature  as  a  function  of  argon  pressure,  together  with  an  estimated  correction. 

The  work  of  Anderson  is  worthy  of  special  note  because  Anderson  compared  probe 
measurements  with  microwave  electron  density  measurements  and  found  wide  discrep¬ 
ancies  between  them.  Anderson  recognized  the  existence  of  serious  perturbation  effects 
in  his  work  and  corrected  the  probe  data  according  to  the  method  of  Davydov  and 
Zmanovskaja.  When  account  was  taken  of  the  influence  of  Coulomb  collisions  on  the 
diffusion  coefficients  the  agreement  between  corrected  probe  values  and  microwave 
values  for  electron  density  was  satisfactory  at  high  electron  densities,  as  shown  in 
Fig.  13.  Also  shown  are  the  ratios  calculated  from  the  Q's  for  various  gas  pres¬ 

sures  with  Coulomb  scattering  neglected.  It  must  be  pointed  out  that,  strictly  speaking, 
neither  the  theory  of  Davydov  and  Zmanovskaja  nor  the  present  theory  ought  to  apply  to 
the  case  studied  by  Anderson,  since  both  assume  a  production  of  ionization  proportional 
to  the  electron  density.  Anderson's  discharge  was  a  glow  discharge  in  which  the  majority 
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Fig.  13.  Data  from  Anderson  on  electron  density  as  measured  by  probes  vs  electron 
density  as  measured  by  microwaves  for  a  cold-cathode  discharge  in  helium 
at  various  pressures.  Also  shown  are  Anderson's  calculations  on  the 
expected  perturbation  according  to  the  theory  of  Davydov  and  Zmanovskaja, 
including  the  effects  of  Coulomb  collisions  ( - )  and  the  expected  pertur¬ 

bation  according  to  these  results,  neglecting  Coulomb  collisions  but  including 
the  fact  that  (O).  For  the  interpretation  of  the  results  at  low  elec¬ 

tron  densities  refer  to  Fig.  14. 


of  the  electrons  were  of  very  low  energy  (T^  ~  500*K)  and  the  ionization  was  produced 
mainly  by  a  very  few  primary  electrons  of  several  hundred  volts  energy  and  by  a  some¬ 
what  larger  number  of  secondary  electrons  with  ~  20,000*K.  A  somewhat  different 
theory  in  which  production  is  assumed  to  be  independent  of  position  would  be  more  appli¬ 
cable  to  Anderson's  case. 

At  low  electron  densities,  Anderson  regarded  his  measured  electron  temperatures 
as  being  too  high,  which  fact  he  ascribed  to  electric  fields  from  the  probe  accelerating 
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electrons  in  toward  the  probe  and  increasing  their  energy.  Too  high  an  electron  tem¬ 
perature  would  lead  to  too  low  a  calculated  electron  density  to  give  the  measured  satu¬ 
ration  current  density. 

Equation  11  indicates,  however,  that  for  negative  sheath  potentials  the  electric  field 

in  the  plasma  is  still  retarding  for  electrons  up  to  zero  sheath  potential.  A  more  likely 

explanation  for  too  high  an  electron  temperature  is  suggested  by  the  type  of  probe  char- 

14 

acteristic  observed  by  Anderson,  shown  in  Fig.  14.  The  presence  of  a  large  electron- 

current  tail,  which  is  due  to  the  presence  of  the  relatively  high-energy  secondary 

electrons,  restricts  the  range  of  currents  over  which  the  temperature  of  the  ultimate 

electrons  can  be  measured  to  the  top  order  of  magnitude  or  so.  Figures  6  and  7  show 

that  the  estimated  electron  temperature  can  easily  be  too  high  by  a  factor  of  from  3  to  10. 

14 

Another  point  to  be  noted  is  that  in  the  internal  report  from  which  his  journal 
article  was  taken,  Anderson  noted  that  the  electron  density  as  measured  from  the 


Fig.  14.  Schematic  of  probe  electron-current  vs  probe -potential  curves 

13 

observed  by  Anderson.  Two-component  electron  energy  dis¬ 
tributions  are  observed:  A,  "ultimate"  electrons  of  very  low 
temperature;  B,  secondary  electrons  (~50,000'K).  Comparison 
with  Fig.  4  shows  that  the  electron  temperature  of  the  ultimate 
electrons  will  be  measured  from  the  range  of  probe  potentials 
near  zero  sheath  potentials,  for  which  the  apparent  electron 
temperature  will  be  too  high. 


32 


electron  current  was  approximately  1 1  times  as  great  as  the  ion  density  measured  from 
the  ion  current.  Since  in  the  plasma  m  =  n^,  the  apparent  difference  must  lie  in  differ¬ 
ences  in  the  degree  of  perturbation.  From  the  differences  between  and  given  in 
Table  I,  we  would  expect  an  apparent  electron  density  approximately  4  times  as  great 
as  the  apparent  ion  density.  It  is  possible  that  the  effect  of  the  shield  could  influence  the 
results  in  the  ion-collecting  range  but  have  smaller  influence  in  the  electron-collecting 
range,  and  thus  further  accentuate  the  difference. 

The  work  of  Bills,  Holt,  and  McClure,  and  that  of  the  author  have  in  common  the 
use  of  a  pulse  technique  for  making  probe  measurements.  Perturbation  effects  will  be 
slightly  different  under  these  conditions,  since  the  pulses  are  too  fast  for  the  density 
distribution  m  the  perturbed  plasma  to  change  during  the  pulse.  The  density  distribution 
in  the  plasma  is  fixed  at  that  determined  by  whatever  bias  is  placed  on  the  probe  between 
pulses.  All  that  changes  during  the  pulse  is  the  electric  field  in  the  plasma.  It  is  easy 
to  show  that  under  these  conditions,  and  under  the  assumption  that  the  probe  is  biased 
strongly  negative  between  pulses,  that 


n^  “  1  +  Q. 
o 


AV  =  - 


Q  (V  +V.) 
V  e  e  1 
''e  “  2(1-H3.)  j 


In  (1+Q.). 


Bills,  Holt,  and  McClure  comment  that  their  electron  density  as  measured  by  pulse 
methods  is  16  per  cent,  and  their  electron  temperature  is  10  per  cent  greater  than  those 
measured  by  dc  methods.  While  it  is  true  that  the  values  of  Q  calculated  for  their  dis¬ 
charge  are  negligibly  small,  their  work  suffers  from  depletion  effects  of  a  different 
sort,  since  the  effective  probe  radius  was  5  per  cent  of  the  tube  radius,  and  the  satu¬ 
ration  electron  currents  were  as  great  as  20  per  cent  of  the  tube  current.  It  seems 
reasonable  to  suppose  that  a  different  sort  of  depletion  theory,  derived  for  this  partic¬ 
ular  case,  would  show  similar  differences  between  perturbations  observed  with  pulse 
and  dc  methods  of  measurement. 
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APPENDIX  I 


I 


MATHEMATICAL  DETAILS  RELATING  TO  SECTION  H 


1 .  Boundary  Condition  on  Electron  Density  at  the  Probe 
From.  Eq.  3, 


n  =  — {a  sin  x  +  B  cos  x} 


X  =  ■ 


irr 


and  from  Eq.  4, 


_1  ^ 

n  dr 


Q 


r 

P 


determine  the  ratio  B/A. 

From  Eq.  5, 

JL  _  JL  2.  ^  -  zf  Acosx-Bsinxj  rt 

n  dr  “  R  n  dx  “  R^Asinx  +BcosxJ  R  x  ‘ 

Equating  this  to  (4),  noting  that  x  «  cos  x  =  1,  sin  x  =  x  ,  we  obtain 

P  P  P  P  ' 


A  -Bx 


Ax  +  B 
P 


>- 


R  x^ 
P  P  P 


A  -  Bx_ 


AXp4B 


1  +  Q 


AXp  -  Bx‘  =  AXp(l-Kl)  +  B(l+Q). 

2 

Since  «  1, 

& 

A=  TTq^p- 

2.  Determination  of  Constant  K  in  Eq.  10 
K 


L 


n  E  +  V  # 
p  p  e  dr 


This  we  evaluate  from  the  probe  boundary  condition  (Eq.  3)  by  solving  for  Ep,  the  elec¬ 
tric  field  in  the  plasma  at  the  sheath  boundary.  Subtracting  the  two  equations  and  sub¬ 


stituting  from  Eq.  4  the  values  and  yields 


n  Q  (V  +V,)  n  Q.(V  +V.) 
p  e  e  i  p  i  e  i' 


=  -2n  E  -  (V  -V,)  ^ 
p  p  e  i  dr 
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Subtract  (V.+V^)  from  both  sides: 
P 


n  Q  (V  +V.)  n  Q.(V  +V.) 
p  e'  e  r  p  i'  e  i' 

- r - + - :: - 


-  (VetV,)  ^1 


-zfn  E  +  V  ^ 
p  p  e  dr 


For  (V  +V.)  If 


^p'^i 

—  (V^+Vi). 

P 


on  the  left-hand  side,  substitute  from  Eq.  4  its  value 


K  =  [nE  +V^ 
I  p  p  e  dr 


n  Q  (V  +V.) 
p  e  e  1 


From  Eq.  8  substitute  the  value  of  n  .  Then 

P 


K 


n  Q  (V  +V.) 
o  e  e  1 

(1+Q)r„ 


3 .  Evaluation  of  the  Integral 


1  +  Q  ’'p  1 


dr 


(l+Q)xJ 


Changing  variables  from  r  to  x  =  gives 

,  r“ 

% 

(l+Q)x 

Now  change  variables  again  to  u  =  — .  Then 

*P 


Qe(V +Vi) 


2 


u 


du 


1  +  Q 


u 


The  integral  can  be  broken  up  by  partial  fractions; 


Qe(V^+V.) 

1+Q  1 


1  +  Q  poo  du 


I. 


1+Q  poo 


u  - 


1  +  Q 


r  ^ 

J,  U 


Q JV  +V  ) 

■■  Q  "in  (1+Q). 
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APPENDIX  II 


MATHEMATICAL  DETAILS  RELATING  TO  SECTION  III 

1.  Derivation  of  the  Ambipolar  Diffusion  Equation  when  —  ^  (rP  )  ~  nv 

r  8r  '  er'  i 

We  assume  that  but,  for  the  moment,  do  not  specify  what  D^  ,  is: 

ar 

-^  =  nv.  +D 

oz  1  ax  r  ar  V  ar /' 


ax 


and  substitute  it  in  Eq.  16. 


^e  -Vj 


(nv.  + 


Ye  -  Y,. 


Yp  -  Y?  D 
e 


1  -  al  )  1  a  /  an\  .  2  B^n 

<V^iV  r  arV^ar^  ^ 


y  +— nv.  =  0 
>ii  1 


>^i<^e"^i)  Ve-\ 

Yf  Yi  *"6 


r  ar  (^a?)  ■*■  '*i<'^e''’^i>  "''i  =  ° 

oz 


The  term  in  front  of  ^  is,  of  course,  D  ,  while  the  term  multiplying  -  ^  (r^)  is  D  . 

az  “  r  or \  ar/  ax 


D  -W^i> 

yf  yf 


Solving  for  D^^,  and  noting  that  »ii(Vg+V.)  is  equal  to  D^,  we  have 
D  = 


p  D 
e  a 


ax  2  2  2 

^i^^e  +  Yef^i 


Multiplying  out  y‘  =  1  +  fxf  B^,  =  1  +  and  noting  that  ^i.  «  gives 


2^2 


D 


D 


D 


1  + 


Therefore 

°a  1  a  /  an\  „  a^n 

„2  r  ar  a?)  +  TT  +  '"'i  =  ® 


which  is  Eq.  20. 
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2.  Solution  of  Eq.  21  Subject  to  Boundary  Conditions  (24) 

_ x/ 

p  dp  8p , 


Assume  that  n  is  a  separable  function  of  position 
n  =  g(p)  f(U 


JL  +  <2  4^2  n 

pg  8p  \P8pj  +  f  a^2  =  0- 

1  I  8^f^  2  2 

Let  -5-1  — =■  =  m.  -  (2,4)  ,  for  which  we  shall  presently  place  a  restriction  on  the  m. . 

^  V847  ^  ^ 

f(4)  =  exp  ^  -  (2.4)^  4^  +  exp^-  -  (2.4)^  4^ 

In  order  for  n  to  be  finite  for  4  =  oo,  all  A.'s  must  be  zero.  Since  g(p)  =  J  (m.p), 

3  3 

n(p,4)  =  ^  ® j  expf-  -  (2.4)^  4)  Jo(m^p)* 

This  must  be  used  as  a  Fourier  type  of  series  to  match  boundary  conditions  at  the  probe. 
Note  that 

P  1 

J  p{Jo(mjp)Jo(mkp)}  »P  *  0  if  j  #  k 


if  j  «  k 


provided  that  m.  and  m.  are  any  two  zeros  of  J  .  Therefore,  we  require  that  the  m.'s 

J  K  O  _  J 

be  zeros  of  J  • 
o 

We  have,  then. 


all  be  zeros  of  j^.  This  also  ensures  that  n  s  0  for  r  s  R  (p=l) 


^  Bj  expf-  -(2.4)^  4)  { 

U*1  J4=0  I 


pUpQ 


pnpQ 


-  Z  ®j  V”?  "5-76Jo(mjP)  = 

i=i 


p  ^  p. 


P  ^  Pr 


P  ^  Pr 


P  ^  Pr 


Multiply  both  sides  by  pj^(m|^p)9p  and  integrate  from  zero  to  unity. 
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Bn  Q 
P  rP 


£  ^  pJo<“kP> 


dp  =  - 


®k  V“k-<2-4)^  jf(nV 


-^-^[mkPpJ^(mj^Pp)]  - - 


nik  -P 


2pn  QJ^(mkP 


®k  =  - 


m, 


-  (2.4)^  jf (nik) 


k  1 


V2  2 

nik  -  (2.4)  =0)  can  be  determined 

by  noting  that  this  term  has  no  ^-dependence,  and  a  p-dependence  of  J^(2.4p).  Its 
coefficient  must  be  the  value  of  n  on  the  axis  an  infinite  distance  from  the  probe,  which 
is  n^,  the  axial  density  in  the  unperturbed  plasma: 

M  exp[-  v/m?  -  (2.4)^  d  Ji(m.p  )  p) 

n(p,U  =  n  J  (2.4p)  -  2pn  Q  Y  - - ; - - 

O  O  P  Lj  /  5  2  o 

j=2  m.  ymj  -  (2.4)^  (m.) 


which  is  the  desired  solution. 


3.  Alternative  Determination  of  Plasma  Potential  Perturbation 
Our  purposes  here  are  to  determine  the  integral 

%Qe<Ve^^i>  r”  dz 
(l+pQ)rp  Jo  n' 

using  a  different  function  for  n,  and  to  show  that  its  form  is  not  greatly  different  from 
the  result  obtained  in  Eq.  31.  Therefore  the  value  of  the  integral  is  not  too  sensitive 
to  the  exact  shape  of  the  function  chosen  for  n. 

It  must  be  pointed  out  that,  since  n  approaches  a  constant  value,  the  integral  will 
diverge.  The  physical  reason  for  this  is  that  the  integral  represents  the  integral  of  the 
electric  field  in  the  plasma  which  is  required  to  draw  the  electron  probe  current  and 
hence  will  ultimately  become  the  integral  of  J^/p^n^.  However,  we  note  that  in  the  pos¬ 
itive  column  of  the  sort  of  discharge  which  we  are  discussing  here,  an  increased  current 
requires  a  higher  electron  density,  but  the  electric  field  remains  approximately  constant, 
or  even  decreases .  In  the  actual  physical  case,  the  drawing  of  probe  current  requires 
that  the  axial  density  at  infinity  be  higher  than  in  the  imperturbed  case.  Since  our  assump¬ 
tions  of  rp  «  R  imply  that  probe  current  will  be  very  much  less  than  tube  current,  we 
shall. neglect  the  disturbance  of  axial  electron  density  at  infinity.  The  divergence  of  the 
integral  will  be  overcome  by  subtracting  the  electric  field  required  to  draw  the  probe 
current  at  constant  axial  density,  n^,  and  write 
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(l+pQ)rp  X 


This  particular  problem  was  avoided  before  by  integrating  between  the  limits  Hp  and  n^ 
this  took  us  not  to  infinity  but  to  z=  prp  with  the  function  chosen. 

PQ  -z/pr^ 


Assuming  that  n  =  -  j  -  -q-  e 


i 


we  obtain 


I  = 


I  = 


"e'  e'  i^  r 

0+PQ)rp  Jo 


PQ  -z/ pr 


-  1 


1 


1  +  PQ 


(1+pQ)  r 


-r 


PQ  J 

1  +  PQ/ 


-z/Pr^ 


PQ  -z/pr 
1 - -  ^ 


dz 


1  +  pQ 

^  PQ^(V^+V.)  p* 

Ju=0  1-Ae-“ 

PQ^+V.) 


ln{l-Ae-X 


1  +  PQ 


PQ  (V  +V  )  f  PQ  I 
1  +  PQ  “  1  +  pQj 
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APPENDIX  III 


CALCULATION  OF  ION  ENERGY  AT  THE  SHEATH  BOUNDARY 


Starting  with  the  fact  that 

V.  -V 
1  g 

Vis  “  —2 - Ep\* 

expressing 


= 


1 


^  “i''c  ”^i  '\/3eV./m. 


3V.m. 
1  1 


and,  from  Eq.  11,  obtaining 


E  =  -V  — 
p  e  n  dr 


Qe(Ve+V.) 


QV^  Q  (V  +V.) 
e  e  e  i 

“'p 

QeVi-SVe 


we  have 


i"  g  /QV  -QV.\/  /3V.m.' 

Replacing  the  Q's  by  their  values,  we  have 
V,  - 


is 


*  —2—  +  V  ^  -P(-VVis) 


V.  /eV 

1  ^  ^ 


V  +V.  p 
e  1  "^e 


Zirm. 


“P<V''es> 


3V.m. 
_ 1 


Here,  the  two  exponential  terms  are  retained  only  when  they  are  less  than  1. 
Multiplying  through,  we  have 

"  'g  V  /3V.  y.  V.  g.  /3V  V.  m. 

''is - -rV7./V  2^.  *m*  '«I><''s''''es>- 

e  1  e  1  e  »  e 
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In  the  vicinity  of  zero  sheath  potential  the  second  term  is  small  in  comparison  to  the 


first;  hence,  for  the  moment,  we  assume  that  «  V^. 


V.  -V 
1  g 

^is  “  2  T 


For  V.  «  V  «  Vg, 


Vis“ 


For  V.  «  V  »  V  , 

leg’ 

V.  V. 

V  = _ i-  +  — i- 

''is  2^2 


n/T 


O  V.  ,  , 

3  ^  1  /n.  J_ 


Therefore,  it  is  reasonable  to  say  that  V.  ,  the  ion  temperature  at  the  sheath  edge, 

xs 

is  approximately  equal  to  the  ion  temperature  in  the  undisturbed  plasma. 
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